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Q_il Abstract 

^^ ■ It is shown that critical phenomena associated with Siegel disk, intrinsic to ID complex 

^vq ' analytical maps, survives in 2D complex invertible dissipative Henon map. Special numerical 

i__i| method of estimation of the Siegel disk scaling center position (for ID maps it corresponds to 

U I extremum) for multi-dimensional invertible maps are developed. 

f^ ■ It is known that complexification of real ID logistic map 

lS; Zn+l = f{Zn) =X~zl, (1) 

^ ■ where X, z E C leads to the origination of the Mandelbrot set at the complex parameter A plane [1] 
r^ and a number of other accompanying phenomena. 

^ I Opportunity for realization of the phenomena, characteristic for the dynamics of complex maps 

■^ . (Mandelbrot and Julia sets etc.) at the physical systems seems to be interesting problem [21 El IH El E] ■ 
In the context of this problem the following question is meaning: Does phenomena of dynamics of 
the ID complex maps (like classic Mandelbrot map ([T])) survive for the more realistic (from the 
point of view of possible physical applications) model - two-dimensional maps invertible in time. For 



X 

;_( ■ example, more realistic model rather than logistic map, is the Henon map 



Zn+l = f{Zn, Wn) = X- Z^-d- Wn, W„+l = g{Zn, ^n) = ^n- (2) 

In the real variable case the system ([21) is 2D invertible map and, hence, can be realized as Poincare 
cross-section of flow system with three dimensional phase space - minimal dimension, providing 
opportunity of nontrivial dynamics and chaos. Henon map is suitable for modelling of the chaotic 
dynamics of the generator with non-inertial nonlinearity, dissipative oscillator and rotator with peri- 
odic impulse internal force etc. [7]. Moreover, Henon map expresses the principal properties of large 
class of differential systems. 

Let us complexify the map ([2D in a such way that z,w, X E C, d eM.. According to the work Pl IS], 
complexified Henon map can be reduced to the two symmetrically coupled real Henon maps and can 
be realized at the physical experiment [4j. 



Let us remark that with \d\ < 1 the Henon map is dissipative system, with \d\ ^ 1 - it is 
area-preserving map, and with d ^ it corresponds to complex ID quadratic map ([T]). 

In previous work [8] we have investigated intrinsic to ID complex maps special scenario of transi- 
tion to chaos through the period multiplication (for example period-tripling) bifurcation cascades and 
found that accumulation points of these bifurcatuins survives for the Henon map with 1 > |(i| 7^ 0. 
In present work we aim to be convinced of existence of one more interesting critical phenomena - so- 
called Siegel disk for the Henon map. This special type of critical dynamics of the one- dimensional 
complex analytical maps corresponds to a stability loss of the fixed point in case of an irrational 
winding number, that is an irrational phase of a complex multiplier with unit modulo /i = e^'^*'^ (for 
example equal to "golden mean" (f = (-\/5 — l)/2). Siegel disk - is a domain inside Julia set at the 
phase plane, filled with the invariant quasiperiodic trajectories rotating around of a neutral fixed 
point. 

Implementation of domain of rotation about the fixed point is connected with existence of local 
variable change z = h{uj) - conjugacy function, which conjugate / with irrational rotation uj —>■ ujh 
near to the fixed point so, that 

fihioo)) = Hfioo). (3) 

This equation, known as Schroder equation, means, that each iteration of / is equivalent to rotation 
of new variable u on an angle 2nip. Conjugacy function (smooth differentiable) exists not at every 
value of complex variable. The domain of its existence corresponds to Siegel disk. The boundary 
of this disk is described by the special fractal quasiperiodic trajectory started from an extremum of 
map (images and pre-images of an extremum are distributed on this boundary dense everywhere). 
Really, let us take derivative of the Schroder equation with regard for h{uj) = z 

f'iz)h'iu)=fih'ifiuj). (4) 

Let's consider a critical value of the variable f'{zc) = 0. From last expression it is evident, that for 
this value (and also for its images and pre-images) conjugacy function should have a singularity - its 
derivatives should turn to zero or infinity.. Besides it is necessary to note, that the extremum is the 
point of a Siegel disk mostly distant from the fixed point. 

At Fig 1 the phase plane for map ([T]) with a value of parameter 

^Siegel ^ 03995409 + 0.5867879^, (5) 

relevant to existence of an irrational neutral fixed point with a "golden mean" winding number is 
represented. Several invariant curves and boundary of a Siegel disk are shown. In a neighbourhood 
of the fixed point invariant curves are close to circles. Invariant curves more distant from the fixed 
point are more distorted. The boundary of a disk, defined by a trajectory, started from an extremum 
is a fractal curve. 



Renormalization group analysis of dynamics at the Siegel disk boundary has been developed by 
M.Widom [S] (see also paper of N.S. Manton and M. Nauenberg [T^). According to this works, two 
types of scaling are possible for quasiperiodic trajectories at Siegel disk: 

1) So-called simple scaling 

which is valid for any point z inside a disk (Here F^- are Fibonacci numbers, and ip - "golden mean" 
winding number); 

2) Boundary of Siegel disk nontrivial scaling law 

f^^+'{z,)-Zc fa, AT -even. 



f^^{zc)-zc { a*, A^-odd, 

with universal scaling factor 

a = -0.22026597 - 0.70848172i. (8) 

Let's consider now complex Henon map. The fixed points of map ((2]) are 

Zl,2 = Wl,2 = • (9) 

There are two complex multipliers /ii^2 for map ([2]), which can be found as eigenvalues of a Jacobi 
matrix 

Jiz,w) = ( ^^(^'")/^^ dgiz,w)/dz \ (" -2^ -d\ ^^^^ 

\ df{z,w)/dw dg{z,w)/dw J \ 1 J 

Thus, multipliers satisfy to the equations 

/ii/i2 = DetJ = d, /ii + Ai2 = TrJ = —2z. (11) 

One can found the value of parameter A, at which one of the fixed points ([9]) is neutral 

A = ^((e'"''^ + rfe-2"'^)2 - 2(1 + d){e^''''^ + rfe"'"'^)). (12) 

Whereas /ii/i2 = d {\d\ < 1), the fixed points with neutral multiplier should have the second multiplier 
with modulo less than unity, i.e. the fixed points should be attractive in directions, orthogonal to 
surface in 4D phase space, at which Siegel disk can occur. 

With d = —0.3 and (f = (-\/5 — l)/2 value of parameter corresponded to a Siegel disk is 

A = 0.0544788851 1456006 + 0.53397699568574151 (13) 

The fixed point is placed in a phase space at 

zq = wo = 0.25807910732741196 + 0.43906869126999041 (14) 

Nontrivial problem is the determination of the scaling center of Siegel disk of Henon map, which 
corresponds to an extremum in the case of one-dimensional map. For this purpose the special 



original numerical method based on universal scaling properties of a Siegel disk ([6]) and ([7]), has been 
developed. 

The main content of a method is following. One should start from any initial point in a multi- 
dimensional phase space close enough to the fixed point (IT^ . After several time iterations of map ([2]) 
trajectory of this point attracts due to a dissipation to one of smooth quasiperiodic invariant curves 
of Siegel disk. Then one should determine the point (z*-^-*, w*-^-*) at this invariant curve, which is 
farthest from rotation center {zo,wo)- One should prolong iterations, until images (of the Fibonacci 
numbers orders) of this point obey with sufficient accuracy to a law of simple scaling. 

According to the nontrivial scaling law of a Siegel disk boundary one can make small step in the 
direction of scaling center and find its approximate location as 

The point {z^'^\ w^"^^) belongs to a new invariant curve, which is farther from the fixed point and closer 

to the boundary of disk. Then one should repeat procedure: 1) find farthest from rotation center 

point of the invariant curve; 2) determine the scale at which this curve is smouth (invariant curves 

more closer to the disk boundary are more distorted), i.e. determine order F^, for which simple 

scaling implements; 3) calculate new approximation of the boundary scaling center using position of 

F/v and -F/v+2 -th iterations of farthest point. 

Repeating procedure more and more times the position of scaling center can be calculated more 

and more precisely {{z^"'\ w^""^) — > {zc, Wc) with n — >• oo). As a result the following position of scaling 

center have been found 

Zc = 0.45756301999 + 0.309658772, 

Wc = 0.22028610776 + 0.64312280?. 

At Fig. 2 the projection of Siegel disk to three-dimensional space is represented. Siegel disk looks 
as smouth surface with fractal boundary. Several invariant curves (concentric lines) on a disk are 
represented. Radial lines at figure correspond to images of a curve, which is a locus of invariant 
curves points mostly distant from the fixed point. Scaling center of Siegel disk is designated. 

The work is performed under partial supoport of RFBR (Grant No 06-02-16619), O.B.I, acknowl- 
edges support from INTAS (Grant 05-109-5262). 
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Figure 1: Phase plane of complex map ([T]). By gray colors the areas of starting values of variable, 
which iterations escape to infinity by different time are marked. Black color corresponds to restricted 
in a phase space dynamics. Curves represent quasiperiodic trajectories around of a neutral fixed point 
(designated by a obelisk). The fractal trajectory started from an extremum of map (designated by 
a circle) represents boundary of a Siegel disk. 




Figure 2: A projection of a Siegel disk for the complexified Henon map ([2]) to a three-dimensional 
phase space. Concentric lines on a disk represent quasiperiodic trajectories. Radial lines correspond 
to images of a curve (bold line with arrow), which is a locus of invariant curves points mostly distant 
from the fixed point. Scaling center of Siegel disk (the point being analog of an extremum) is 
designated as {zc,Wc). 



